Measurement of the Coulomb Logarithm in a Radio-Frequency Paul Trap 
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Samples of ultracold 174 Yb + ions, confined in a linear radio-frequency Paul trap, are heated via 
self-induced micromotion interruption, while their temperature, density, and therefore structural 
phase are monitored and simulated. The observed time evolution of the ion temperature is compared 
to a theoretical model for ion-ion heating allowing a direct measurement of the Coulomb logarithm 
in a linear Paul trap. This result permits a simple, yet accurate, analytical description of ion cloud 
thermodynamic properties, e.g. density, temperature, and structural phase, as well as suggests 
limits to and improvements for on-going trapped-ion quantum information efforts. 
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Binary collisions in an ensemble of charged particles 
are fundamental throughout physics. As such, modeling 
their behavior plays an important role in applications 
ranging from thermonuclear fusion to quantum compu- 
tation. Despite their importance and the large amount 
of work towards understanding their effects, there is still 
considerable ambiguity in how to best model these colli- 
sions. Since Landau's early work [![, the most straight- 
forward approach handles the divergence associated with 
collisions of charged particles by introducing both a 
short-range and long-range cut-off for the 1/r interaction 
potential. The long-range cut-off is typically associated 
with the Debye screening length, Ad = \J e ksT / pe 2 , 
while the short-range cut-off is associated with the Lan- 
dau length, Rc = e 2 /(4ire kBT), which is the distance of 
closest approach. Taken together, the integration bounds 
give rise to the so-called Coulomb logarithm in its sim- 
plest form, lnA = In (CX D /R C )- 

Over the last 75 years, there have been many at- 
tempts to calculate an accurate form of In A, ranging 
from straightforward estimations of the coefficient C [2| 
to sophisticated analytical treatments Q with reasonable 
consensus that C « 0.765 These results give satis- 
factory agreement with data for weakly coupled systems, 
g = Rc/^d <C 1, but clearly fail to describe strongly- 
coupled systems, g ^> 1, where the collision rate satu- 
rates. In this regime, more sophisticated treatments Jj|, 
which do not require a short-range cut-off, have recom- 
mended several alternative forms for In A. And recently, 
a new approach, motivated by the need to model ther- 
monuclear ignition, used molecular dynamics simulations 
to suggest In A - In (1 + 0.7 /g) @, for g < 10. 

Given the importance of the Coulomb logarithm, there 
have also been attempts H, H at a direct measurement 
of its dependence on the strong-coupling parameter g. 
However, to date, these experiments have either been 
confined to the weak-coupling limit or were inconclusive 
because of the difficulties associated with parameterizing 
dense, high-energy plasmas. 

In this work, we use an alternative route to realize 
a strongly-coupled ion system and measure In A: laser- 



cooled 174 Yb + ions confined in a linear Paul trap. Here, 
the confining trap potential provides a smoothly vary- 
ing, neutralizing background for the positively charged 
ions, resulting in a system described as a one-component 
plasma. Though the density of this plasma is orders of 
magnitude lower than a typical strongly-coupled plasma, 
the low temperatures accessible through laser cooling 
make it possible to realize systems with g 1. Fur- 
ther, by laser cooling the sample to a large g and then 
allowing the ions to heat through self-induced micromo- 
tion interruption, we are able to measure the evolution of 
both the temperature and structural phase of the trapped 
ion cloud over a large range in g. From these measure- 
ments, we are then able to determine the Coulomb log- 
arithm over a range of 10" 7 < g < 10~ 2 . Using a de- 
tailed molecular dynamics simulations, we confirm this 
experimental determination and extend it to a range of 
10 -7 < g < 10 3 . As the values of In A for large g are 
known to be process dependent (l(]| . this result must be 
carefully interpreted before it can be applied to other sys- 
tems. Nonetheless, it offers a complete description of In A 
for radio-frequency Paul traps, and thus allows a simple 
analytical description of trapped ion thermodynamics. 

In the remainder of this manuscript, we explain the 
phenomenon of ion heating by self-induced micromo- 
tion interruption, detail the method by which In A is 
extracted, and describe the experimental system. We 
present experimental and molecular dynamics results and 
a recommended expression for In A in linear Paul traps. 
We conclude with a discussion of the implications of this 
work for trapped-ion quantum information efforts. 

In a linear Paul trap, an ion experiences both a time- 
dependent force from the confining electric potential of 
the trap and Coulomb repulsion from the other ions, re- 
sulting in trajectories given by: 
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with the trap potential given as: 
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where r Q is the field radius, z is the distance from trap 
center to the end cap used for axial confinement, V r f and 
V ec are the radio-frequency (rf) and end cap voltages, 
respectively, fl is the frequency of the rf voltage, and a 
is a geometric factor less than unity. Due to the infinite 
range of the Coulomb interaction, Eq. ([1]) represents 
a complicated many-body problem and has no closed- 
form solution. Therefore, two alternative approaches are 
usually employed to treat this system. 

First, rigorous molecular dynamics (MD) simulations 
have been performed to study the structure of ion clouds 
[IH and rf heating rates and for comparison with 
experimental ion fluorescence images Despite their 
successes, these simulations offer little physical intuition, 
making it difficult to optimize a given system. 

The second approach to modeling large ion systems 
in a Paul trap has been through analytical techniques 
l44l7l|. Of the analytical approaches, the simplest and 
most intuitive describes the trapped ion trajectories by 
the well-known Mathieu solutions and includes the ef- 
fect of the Coulomb interaction as hard-sphere colli- 
sions between the trapped ions to calculate, among other 
things, the evolution of the trapped ion kinetic energy 
16L 17 1 . In this limit, the kinetic energy evolution of an 
ion due to the collisions with other ions is given simply 
as, W — jAW, where 7 is the collision rate and AW is 
found from the kinetic energy change per collision [l?} af- 
ter enforcing conservation of momentum and energy for 
the collision and requiring that the new ion trajectory 
corresponds to a Mathieu solution. In contrast to static 
traps, AW does not average to zero over time, or over 
the ensemble, in a rf Paul trap. In fact, as shown in Ref. 
[13] upon averaging AW is always positive, leading to the 
so-called micromotion-intcrruption heating phenomenon. 
Though this heating has been explained in different ways 
3 3 111 , it arises from the simple fact that when ions 
undergo collisions their trajectories are not given by the 
stable Mathieu trajectories and as a result the rf trapping 
field can do net work on them. 

To link temperature with W , we introduce T sec and 
T to t to be proportional to the random thermal energy 
(secular motion) of the ions and the total kinetic energy 
(secular motion plus micromotion), respectively [l9| : 



CT U 



W 



CT sec = W sec = V (W - W ex ) (3) 



where C = |7Vfcs, 77 is the ratio of secular energy to total 
energy and ?7 ~ § for low Mathieu q parameter as a con- 
sequence of equipartition of energy between secular mo- 
tion and micromotion |l6l|. W ex accounts for the excess 
micromotion energy [20( due to displacement of the ion 
from the node of the oscillating electric field, as a result 



of cither the location of the ion in the crystal or stray, un- 
compensated, dc electric fields. Typically, T ex = W ex /C 
is a few Kelvin for an ion crystal composed of N ~ 10 3 
" H. 
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Using the Chandrcshakar-Spitzer plasma self-collision 
rate Q , the rate of change of the total and secular tem- 
perature of the ion cloud is given as: 



Tt, 



e 4 pj{T sec ) In A 
47re§ v ^(fc JJ T sec )3/ 2 
e A pi(T sec ) In A 



tl,, 



e{T sec + r]T ex ) (4) 



where Pi(T sec ) is the ion density [13] and e = AW/W is 
the average fractional increase of the ion energy per colli- 
sion. By averaging over the rf-phase at which the collision 



takes place, Refs. [la, [17J have calculated e in terms of 
the Mathieu stability parameters a and q. Through nu- 
merical integration of their result, we have found e can be 
simplified to e « |(l+2g 2 - 24 ) with a relative error < 0.4% 
for q < 0.4, a = 0. Thus, by laser cooling a sample of 
trapped ions to a low initial temperature, extinguishing 
the laser cooling, and monitoring the ion temperature 
evolution, we are able to measure In A as a function of g. 

The experimental apparatus used in this work consists 
of a sample of 174 Yb + ions, loaded via laser ablation, into 
a linear rf Paul trap with r a = 1.2 cm, z a = 1.075 cm, 
Tj = 0.13, n = 2tt x 300 kHz, V rf = 155 V and V ec = 
5 V. A strongly coupled ion ensemble (N = 10 2-3 ) is 
realized by laser-cooling the Yb + ions, along the trap 
axis, with a 369 nm cooling laser (detuned from res- 
onance by 8 = —30 MHz) and 935 nm repump laser 
(5 = MHz) to a starting secular temperature, mea- 
sured from the Dopplcr broadened fluorescence profile, 
ranging from Doppler limited, Td ~ 1 mK to 100 mK, 
depending on crystal size, resulting in a one-component 
plasma with g = 10 2 ~ 10 3 . 

Once the strongly-coupled plasma is established in the 
trap, the cooling laser is extinguished and the ions evolve 
in the trap and heat through self-induced micromotion 
heating. After a variable time delay, the cooling and re- 
pump lasers are reapplied and the fluorescence level of the 
ion cloud immediately recorded. If the ion temperature 
has increased during the time when the lasers were extin- 
guished, this fluorescence level will be different than the 
steady-state value reached for the initially cold plasma, 
sec Fig. QIa). By recording the ratio of fluorescence be- 
fore and after heating, the temperature of the ions can be 
estimated as a function of heating time in a manner sim- 
ilar to Ref. 2l|, as shown in Fig. []Jb)-(c) for a sample of 
ions with N = 280. Typically, the observed fluorescence 
ratio decreases with increasing temperature since both 
the fluorescence profile is further Doppler broadened and 
the higher energy ion trajectories have less overlap with 
the laser beam inducing the fluorescence, see Fig. Ula)- 
(b). However, for finite laser detuning the fluorescence 
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FIG. 1: (a) Laser fluorescence profile for a sample of ions at 
the Td (solid) and at ~90 K (dashed). The arrow denotes the 
change in fluorescence at a fixed detuning of 8 — -30 MHz. (b) 
The observed (dots) and simulated (line) fluorescence ratio for 
S — -30 MHz vs. heating time, (c) The extracted (dots) and 
simulated T sec (line) and T to t vs. heating time. 



ratio may actually increase for a very small temperature 
range just above the Doppler temperature. In this range, 
the fluorescence profile is slightly broadened so that the 
detuned laser addresses more ions. This effect gives rise 
to the slight peak in fluorescence ratio seen at .03 s in 



Fig. QJb); similar features were observed in Ref. 21 1. 

As seen in Figs. QJb)-(c), the fluorescence level, and 
thus the temperature, is relatively unchanged until 10- 
100 ms after the laser cooing is extinguished, when a 
sharp increase in temperature occurs, followed by a re- 
gion of slower heating. As detailed below, the relatively 
small heating rate observed at early times is a conse- 
quence of the suppression of ion-ion collisions, i.e. a small 
In A, for a strongly-correlated plasma, while the sudden 
jump in temperature coincides with the phase transition 
from the liquid to gas phase. At the liquid-gas boundary, 
the ion density is still relatively high, but the ion motion 
becomes less correlated, i.e. increased In A, leading to a 
larger heating rate. As the ions move into the gas phase, 
the motion becomes even more uncorrelated, leading to 
a further increase in In A, however, the density, and thus 
the collision rate drops, leading to a reduced heating rate. 
Also shown in Figs. HJb)-(c) are the results of a molecu- 
lar dynamics simulation, which initializes the ions at the 
experimentally realized temperature and then integrates 
Eq. ([l} numerically using a leap frog algorithm [22j im- 
plemented in ProtoMol software [23j. As the ions heat 
through self-induced micromotion interruption, their flu- 
orescence level is calculated from the known laser inten- 
sity profiles and a rate-equation model, which includes 
the variation of laser intensity and Doppler shift for each 
ion position and velocity, repectively. Given experimen- 
tal imperfections, such as stray fields, machining errors, 
laser amplitude and frequency noise, etc., that are not 
included in the simulation, the agreement between the 
simulated and measured fluorescence ratios (Fig. [lib)) 
is satisfactory. In what follows, we use these results to 
extract the ion-ion heating rate and ultimately In A. 

Using experimental data like that shown in Fig. [1] for 
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FIG. 2: The experimental (black dots) and molecular dy- 
namics (white dots) determinations of In A versus g. Despite 
large variation in trap parameters (see text) the observed val- 
ues fall along the same curve, indicating a 'universal' form for 
In A. The red line represents the best fit described in the text, 
while the black and dashed line are the results of Ref. [(| and 
ln(0.765/<?) [H, respectively. Shown at the bottom of the fig- 
ure are the pair correlation functions and velocity distribution 
functions, extracted from simulation, for selected p-values. 



ion clouds with N between 300 and 3000 and modeling 
the density as 



Pi (^sec) 



T < T 
-l sec _^ p 




37V mr^tt 
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where p max = 4e K/ / mr^n 2 , ± p - 2ki 
and U is the geometric mean of the three secular frequen- 
cies, Eq. [4] is inverted to find In A and the results are 
plotted in Fig. 2. As the heating rate is the time deriva- 
tive of T sec , the coarse granularity of the experimental 
data in time makes it difficult to calculate reliable values 
of In A at short time scales. Therefore, we also deter- 
mine the heating rate and In A from molecular dynamics 
simulations. For the simulated data, the heating rate 
is found by taking the numerical time derivative of the 
ion temperature, defined by \NkBT to t = y Y^iLi ^i(*) 2 
and \Nk B T sec = y X)»=i v sec ,i(t) 2 where Vi and v aec ,i 
are the total and secular velocity of i th ion and the over- 
line denotes averaging in principle over t = {—00,00} - 
in practice averaging over several secular motion periods 
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is sufficient. The resulting values for In A are consistent 
with those extracted from experimental data, as shown 
in Fig. [3J but are expected to be of higher accuracy. Us- 
ing this technique, molecular dynamics simulations were 
performed, like those shown in Fig. [TJ for a range of 
ion numbers TV = {50, 100, 500, 900} and ion-cloud ge- 
ometries (radial-to-axial cloud size aspect ratios 24| of 
R/z = {0.25, 1,4}) to determine if the parameterization 
of Eq. E] leads to a universal form for In A in Paul traps. 

The values of In A extracted from the simulation are 
plotted versus g in Fig. [2] alongside the Landau-Spitzcr 
result [3] and the result of Ref. 0] . Also shown, as the 
top x-axis in this figure, is the corresponding plasma cou- 
pling parameter Y = e 2 / (Aire a ws kBT sec ), which, given 
the Wigner-Seitz radius a ws = {/?>/ {^npi) and secular 
temperature, characterizes the structural phase of the ion 
cloud as denoted by the three regions of the graph [25| . 
Clearly, despite the large changes in ion number and ion- 
cloud geometry, the dependence of In A on g appears uni- 
versal and can be parametrized by the piecewise fit: 



In A 



ln(l+0.7/g) 
1+500 

fu(g) = ^ 



S<1 
5>1 



(0) 



where the form of fi(g) has been inspired by Ref. [6|. 
Interestingly, the observed change in dependence of In A 
on g occurs near the gas-to-liquid phase boundary, which, 
since Xd/cl ws = l/\/3<?, also coincides with the regime 
where the Debye length becomes smaller than the average 
inter-particle spacing. Therefore, assuming Debye theory 
is approximately valid for g > 1/3, it is reasonable to 
expect the ion-ion cross-section is proportional to Xjj and 
thus In A cx XjjT 2 oc g~ 2 in agreement with the fit. 

For reference, the structural phases of the ion cloud 
are presented in Fig. [2] by the pair correlation function 
[IK G(r/a ws ) = (Anp l r 2 y 1 dN/dr for three different g 
values, as determined by molecular dynamics simulation. 
For the solid phase G(r/a ws ) exhibits a sharp peak at 
r/a ws ~ 1.7, confirming a highly-ordered crystal struc- 
ture, which disappears as the ion cloud moves into the 
liquid and gas phases. Also, shown at the bottom of 
Fig. [5] are the secular velocity distributions of the ions 
at the same three g values. In this figure the points are a 
histogram of the simulated secular velocity distribution, 
while the solid curve is the Maxwell-Boltzmann (MB) 
velocity distribution expected for the calculated temper- 
ature. Clearly, the simulation results for the gas and liq- 
uid phases are consistent with the MB distribution, con- 
firming the appropriateness of the Chandreshakar-Spitzcr 
rate in deriving Eq. @] For the solid phase, the velocity 
distribution exhibits a significant power-law tail, violat- 
ing the assumptions of Eq. [4] and preventing an accurate 
determination of In A in this phase. 

Several analytical results can be derived using Eq. [6] 
that provide insight into the plasma dynamics and have 
important consequences for quantum computation with 
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FIG. 3: Comparison of experimental temperature, density 
and heating rate (dots) for two ion clouds (iV = 280 (blue) and 
N = 2800 (red)) as a function of time to the result predicted 
by Eq. [3] with In A given by Eq. [S](line). 



trapped ions. First, for g <C 1, In A can be approximated 
as a constant (as is typical for low density plasma [26l |) 
and Eq. ^integrated, yielding T sec oc i 1 / 3 . In this regime, 
the ion temperature grows slowly until it is eventually 
balanced by evaporative loss from the trap or sympa- 
thetic cooling from residual neutral background gas. Sec- 
ond, assuming that T ex T sec in the initial ion crystal, 
Eq. [4] can be directly integrated for g > 1 yielding T sec = 
[c(i m -t)]~ 2 , where c = fi(l)2Tre k B /2 er]T ex /(y/rne 2 ) and 
the time of the dramatic rise in temperature as the ions 



move into the gas phase is t T 



-1/2 



(t = 0)/c. As 



this expression for t m ignores the contribution of T sec to 
the heating rate, it slightly overestimates the time of the 
phase transition; however, in the case of our experiment 
we have found it accurate to within ~ 10 ms. Interest- 
ingly, in quantum computation with strings of trapped 
ions the computational gate operations occur with the 
laser cooling extinguished. Therefore, t m represents the 
upper limit for the time to implement a computational 
algorithm since, once the ion string melts, the quantum 
information is lost and the register must be reinitialized. 
For the parameters of e.g. Ref. j2?| with 14 ions, we 
find this fundamental limit to be ~ 10 3 s. If this sys- 
tem is scaled to a larger number of ions, as necessary for 
many practical quantum computation applications, t m 
will be signficantly reduced if excess micromotion is not 
controlled and may limit the number of possible gate op- 
erations. Likewise, recent proposed experiments to use 
kinked ion chains to study the Kibble-Zurek mechanism 
[28| and the coherence of discrete solitons [29| will be 
fundamentally limited to timcscales less than 1 ~ 10 s. 
In addition to providing the upper limit for a single com- 
putation/simulation, the expression for t m can be used 
to guide future efforts. For example, linear string geome- 
tries of heavy ions at low Mathieu q parameter should 
exhibit the longest lifetimes. 

Finally, to demonstrate the utility of the expression 
for In A, Fig. [3] compares experimental data for two ion 
clouds of different size (N = 280 and 2800), taken in the 
same manner as the data of Fig. [TJ with the temperature 
and density predicted by the integration of Eq. U using 
Eq. [5] In addition to providing a simple means to ac- 
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curatcly calculate the thcrmodynamical properties of a 
system of trapped ions, these expressions explain several 
well-known experimental observations. For example, the 
smaller ion-ion heating rate (Fig. (He)) in the solid phase 
is due to ion-ion correlation as quantified through In A. 

In conclusion, we have measured the heating rate of 
ions trapped in a linear rf Paul trap due to self-induced 
micromotion interruption. These data, and detailed 
molecular dynamics simulations, have been used to de- 
termine the value of the Coulomb logarithm over a range 
of 10" 7 < g < 10 3 . Though most determinations of In A 
are process dependent [10j, we expect our results to be 
comparable to In A in other one-component plasmas for 
g <C 1, as the Mathieu trajectories accurately describe 
the ion motion in this regime. This expectation is sup- 
ported by the fact that our result converges to the tradi- 
tional Landau-Spitzer result in this regime. However, as 
g grows the Mathieu solutions provide a less accurate de- 
scription of the ion trajectories, leading to a change in e.g. 
e. If future theoretical work accounts for these effects, 
then our measurement might be reinterpreted to give a 
model independent determination of In A. Nonetheless, 
in its current form our result permits a simple, yet ac- 
curate, analytical description of ion cloud temperature, 
density, and structral phase transitions in a linear Paul 
trap. Thus, it should be immediately useful to a number 
of experimental efforts, including the growth of large ion 
crystals 13011 , sympathetic cooling of atomic or molecular 
ions [3ll433l | , and trapped- ion quantum information. 

This work was supported by ARO grant No. W911NF- 
10-1-0505 and NSF grant No. PHY-1005453. 
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